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Abstract: Recently, the D0 Collaboration measured the CP- violating like-sign 
dimuon charge asymmetry in neutral B decays, finding a 3.2a difference from the 
standard-model (SM) prediction. A non-SM charge asymmetry a s sl suggests a new- 
physics (NP) contribution to B s -B s mixing. In this case, in order to explain the 
measured value of a s sX within its la range, NP must be present in rf 2 , the absorptive 
part of the mixing. In this paper, we examine whether such an explanation is possible 
in models with flavor-changing Z (ZFCNC) or Z' (Z'FCNC) gauge bosons. The 
models must also reproduce the measured values of the indirect CP asymmetry 
in B s ->• Jfij)<j>, and Ar s , the B s -B s width difference. We find that the ZFCNC 
model cannot reproduce the present measured values of and within their la 
ranges. On the other hand, in the Z'FCNC model, the values of all three observables 
can be simultaneously reproduced. 



KEYWORDS: Flavor- Changing Gauge Boson, B s -B s mixing, Beyond Standard 
Model. 
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1. Introduction 

Flavor-changing neutral-current (FCNC) processes occur only at the loop level in the 
standard model (SM), and are therefore a very sensitive probe of physics beyond the 
SM. Interestingly, recently there have been several measurements of quantities in B 
decays, especially in the b — > s transition, which differ from the predictions of the 
SM. For example, i) the value of sin 2/3 measured in several b — > s penguin decays 
is found to be smaller than that obtained in Bj — > J/ipKs [0, ii) in B — > ttK, the 
SM has some difficulty in accounting for all the experimental measurements |2|, iii) 
the forward-backward asymmetry in B — > K*[i + fi~ shows a small deviation from the 
SM, especially at low momentum transfer squared [JJ], iv) the measurement of CP 
violation in B s — > J/ip^> suggests the presence of a new-physics (NP) CP- violating 
phase in B s -B s mixing || . Although these effects are not significant enough to claim 
NP, at least they indicate that flavor physics will still play a very important role in 
the search for NP in the LHC era. The LHCb will provide much more precise data 
on the above-mentioned observables. 

Recently the D0 Collaboration reported an anomalously large CP- violating like- 
sign dimuon charge asymmetry in the B system 0. The asymmetry is 

^ s "£-*r , (i.i) 
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where iV^ =± is the number of events of bb — > fj^fj^X. From the analysis of 6.1 fb 1 
of data, D0 found 

A b sl = -(9.57 ± 2.51 ± 1.46) x 10~ 3 , (1.2) 

where the first (second) error is statistical (systematic). This exhibits a 3.2a devia- 
tion from the SM prediction, A b s [ SM = (— 2.3^'^) x 10~ 4 0. The asymmetry can be 
written as a linear combination of and B s contributions 

A b A = (0.506 ± 0.043)4 + (0.494 ± 0.043)4 , (1.3) 

where the semileptonic "wrong-charge" asymmetry o'j (q = d, s) is given by 

aQ = T(B q -» n+X) - Y{B q -» i£X) 
sl T(B q -> + -)■ ' 

Thus, the discrepancy may indicate a NP contribution to B^-B^ mixing and/or B s -B s 
mixing |7|]. 

Using the current experimental value of a^ exp = — 0.0047 ± 0.0046 || |8| and the 
SM prediction a d f M = (-4.8lJ;§) x 10" 4 , Eq. Q leads to a", = -0.0146 ± 0.0075. 



The CDF Collaboration also measured -A^, but with much larger errors: = 
(8.0 ± 9.0 ± 6.8) x 10~ 3 §. Finally, D0 directly measured a s sl : af xp = -(1.7 ± 
9.ll^; 4 ) x 10~ 3 110]. Combining all these results, we get 

K s i)avo = -(12.7 ± 5.0) xlO- 3 , (1.5) 

which is still about 2.5a away from the SM prediction of a'f = (2.1 ± 0.6) x 10" 5 . 

The D0 measurement thus suggests the presence of NP in B s -B s mixing, and 
we explore this possibility below. We begin with a general review of the mixing. In 
the B s system, the mass eigenstates Bl and Bu (L and H indicate the light and 
heavy states, respectively) are admixtures of the flavor eigenstates B s and B s : 

\B L ) =p\B s )+q\B s ) , 

\B H ) =p\B s )-q\B s ) , (1.6) 

with |p| 2 + |q , | 2 = l. As a result, the initial flavor eigenstates oscillate into one another 
according to the Schrodinger equation 

4(\m)=(M--i r 4)(\m) , (i.7) 



dt\\B,(t))J \ 2 J \\B,(t)) 

where M = and r = correspond respectively to the dispersive and absorptive 
parts of the mass matrix. The off-diagonal elements, M[ 2 = M|* and T\ 2 = F^*, are 
generated by B s -B s mixing. We define 

r s = TH t T \ AM S = M H -M L , AT S = T L -T H . (1.8) 
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Expanding the mass eigenstates and q/p in Tl 2 /M( 2 , we find, to a very good approx- 
imation fill, 



AM S = 2\M* 2 \ , 

Ar s = 2|r^ 2 | cos0 s , 



1- ^-i4sin0 s , (1.9) 

v V 2|r s 12 | 

where 0^ = argMf 2 and 4> s = arg(— Mf 2 /r^ 2 ). Then is given by 

a- =Imf^ =£M s in0 s . (1.10) 



In the SM, we have 

a s sl = (4.97 ±0.94) x 10~ 3 sin0 s . (1.11) 

Within the SM, the phase of Mf 2 is given by arg[(V#,VJ*) 2 ] and that of rf 2 is domi- 
nated by arg[(V cb V; s ) 2 ]. The state-of-the-art calculation gives <p s = 0.0041±0.0008 [|. 
As noted above, the SM cannot explain the D0 result; this is due to its far-too-small 
weak phases. From here on, we neglect the SM weak phases. 

Suppose now that NP contributes to M[ 2 but not to rf 2 . We introduce a new 
parameter A s to take this into account: M( 2 = M* 2 SM + M^ 2 NP = M^ 2 SM A S = 
M° 2 S M |A.Je^M. We have 



< = T§m s -ix¥ = < 4 - 97 ± °- 94 ) x 10 " n^f • (I-") 



The experimental value of the mass difference in the B s -B s system, AM S = (17.77 ± 



0.12) ps _1 , constrains |A S | = 0.92 ± 0.32. Then Eqs. ([OjD and ( |rT2|) give sin 



—2.56 ± 1.16, implying that the full la range of the experimental result lies outside 
the physical region. It is therefore not possible to explain the la range of the D0 
measurement if NP contributes only to M( 2 . 

This problem can be solved if the new physics contributes to T\ 2 . Here we 
introduce a second new parameter E s : T\ 2 = r^ 2 SM + r^ 2 NP = rJ 2 SM S s = r^ 2 SM |S s |e^r. 
We then have 



ipa.SMi 
I 1 12 






|M x f M 




A s | 



«sl = llrt .SM,IA i sm ^M 



= (4.97 ± 0.94) x 10~ 3 i^L sin(0^ - 0f ) . (1.13) 

I s I 

Now Eq. ( p~~5|) can be reproduced if |5 S | is sufficiently large. 

The difficulty is that, in many NP models, the contribution to Tf 2 is not large 
enough to compete with that of the SM, which is dominated by the tree-level b — > 
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sec process. Two exceptions are R-parity- violating supersymmetric models [12] and 
leptoquark models [|13| (other analyses of Eq. (|1.5| ) in various NP models can be 



found in Refs. 0, III). 



In this paper we examine the contribution of flavor-changing neutral gauge 
bosons to T s l2 . We consider two types of models. They involve tree-level Zbs or 
Z'bs couplings. As we will see, the non-universal Z' model can enhance enough 
to explain the la range of Eq. ( |1.5p . 

The paper is organised as follows. The NP contribution to TJ 2 is discussed in 
Sec. 0. In Sec. [3], we describe the models with flavor-changing Z and Z' couplings. 
The numerical results are presented in Sec. |]. We conclude in Sec. |5|. 



2. Contributions to 



In the SM, the dominant contribution to r^ 2 SM comes from the charged- current b — > 



sec operator with a c-quark loop. It is given by 

■ns,SM 



12 



3tt 



K 1 (l-x c ) + -K 2 {l-Ax c 



(0 LL ) + {K t - K 2 )(l + 2x c )(O rr ) }, (2.1) 



where A r 



VcbV* x c 



ml /ml, and K\ 



3C? + 2dC 2 , K 2 



cl 



The values for 



the Q are d(m 6 ) = 1.086, C 2 (m b ) = -0.197 (Tables 4, 5 of Ref. @). We use the 



vacuum insertion approximation to calculate the hadronic matrix elements [16 
(Oll) 

(Orr) = (B s \sP R b sP R b\B s ) = ~m B J% t ( J^l ) , (2.2) 



[B s \s ltl P L bs^P L b\B s ) = l -m Bs f Ba 



mb + m s 



where P LjR = (1 =F j 5 )/2 and f Ba = (238.8 ± 9.5) MeV 0. An updated theoretical 
prediction can be found in Ref. , whose results we use in our numerical calculations. 

Any NP contribution to must come from a new operator of the form b — > 

sf f, where / is a light fermion. If r^2 NP is to be significant, it must be at least 
comparable to r^ 2 SM . Now, most NP light fermionic operators are constrained to 
be small. In particular, (i) e and \i loop contributions are strongly bounded by the 
b — > se + e~ and b — >■ s[i + [i~ processes, (ii) the b — >■ suu and b —> sdd operators 
are constrained by the measurement of B[B — > ttK), (iii) B{B — > <pK) constrains 
b — > sss. On the other hand, the bounds on the NP b — > sr + r _ transition are very 
weak. For example, the present upper bound on B(B S — > r + r _ ) is only < 5% |nj, 
to be compared with the SM prediction of B(B S — > t + t~) ~ 10 -8 . Also, the current 
upper bound on B{B — \ X s t + t~] 
models we consider, we examine the contribution to 



is just < 5% [0. For this reason, in the NP 

coming from b — > st + t~ . 
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Another potential important effect on T s {^ p comes from the NP b — > sec operator. 
Although the constraints from B — > DD S , Bd — > J/ipKg, etc. are such that the NP 
c-quark loop contribution cannot be large enough to compete with T^ SM , the SM-NP 
interference term can be significant. 

The decay width difference Ar s has been measured independently. The angular 



analysis of B s ->■ J/i(j(j) gives f8|, |I9|, [20 



Ar s = ±(0.154±g5§*) ps- 1 , (2.3) 

to be compared with the the SM prediction |J 

AT™ = (0.096 ± 0.039) ps" 1 . (2.4) 

If NP contributes to rj 2 NP , it is present in the width difference, whose expression is 
given by 

Ar s = Arf*|H s |cos(<^-0f). (25) 

The measurement of Ar s therefore constrains Fi 2 NP . 

There is another measurement which must be taken into account. The CDF [Ell 



and D0 Collaborations have measured indirect CP violation in B s — > J/ipcfr. 
They obtain = —2(3 S , and find || 

(5 S = 0.39i°i« or 1.18^ . (2.6) 

This disagrees with the SM prediction 

/3 S SM = 0.019 ± 0.001 (2.7) 

at 2a. (Note: the recent CDF measurement agrees with the SM better than the 
previous D0and CDF measurements.) Now, we have assumed that the NP affects 
M( 2 , so that there is a weak phase in B s -B s mixing, <p s M . Depending on which NP 
operator(s) contribute to r^2 NP , there may also be a contribution to the decay of 
B s — > J/ip<f). Thus, the NP is constrained in that the CDF/D0 measurement must 
be reproduced. 

To detail the contribution of the NP to the indirect CP asymmetry in B s — > 
J/ip4>, we follow the procedure of Ref. [p3]| . First, we note that B s — > J/ip<p is really 



three separate decays, one for each polarization state A of the final-state vector 
particles; longitudinal: A = 0, transverse: A = {||,_L}. Second, in Ref. |^3| it is 



argued that all strong phases associated with NP amplitudes are negligible. In this 
case, for each polarization one can combine all NP matrix elements into a single NP 
amplitude, with a single weak phase (p\: 



((W)a| O np \B s ) = b x e^ . (2.8) 
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We now assume that this single NP amplitude contributes to the decay of B s — > 
J/ip(j). The decay amplitude for each of the three possible polarization states may 
then be written as 

A\ = Amp{B s -> J/^) A = axe^-W + bxe^e'^ , 

Ax = Amp{B s ->■ J/ip<j))x = a A e <( *X-*D + bxe'^e'^ , (2.9) 

where a a and 6a represent the SM and NP amplitudes, respectively, ipx is the new- 
physics weak phase, and the 5" are the SM strong phases. All strong phases are given 
relative to 5j_. ax is defined to be positive for every polarization. 6a can also be taken 
to be positive: if it is negative, the minus sign can be absorbed in the weak phase 
by redefining (fx — > tpx + it. We emphasize this fact by writing the ratio bx/ax as the 
positive-definite quantity \r\\. Note that strong phases are generated by rescattering, 
and this costs a factor of about 25. The strong phase of the SM color-suppressed 
b —> ccs diagram C is generated by rescattering of the color-allowed b — > ccs tree 
diagram T. Since \C/T\ is expected to be in the range 0.2-0.6, the SM strong phase 
is on the small side, but is not negligible. 

Putting all this together, we find that the indirect CP asymmetry in B s — > J/i/xf) 
measures 

= sin <f>* M -\-2\rx \ cos (f> s M sin <px cos 5" • (2.10) 

(We have neglected the SM contribution to the weak phase of B s -B s mixing, which 
is expected to be only ~ 2% [Eq. (|2.7|) 1.) If the NP contribution to comes 
only from b —> st + t~, then r A = and = smcj) s M . However, if the NP b — > sec 
operator is involved in Ti^, then rx ^ and the full expression for above must 
be used. 

Finally, we note that the phases of a|j and are different when there is a NP 
contribution to r^ 2 and/or B s — > J/ipcf). In addition, the relation between them [ f25|j . 

IAr, 



I AM, 



is violated in this case. 



3. New-Physics Models 
3.1 Z-mediated FCNC's 

In the model with Z-mediated FCNC's (ZFCNC), a new vector-like isosinglet down- 
type quark d! is added to the particle spectrum [^6|, ^7], Such quarks appear in 
Eq GUT theories, for example. The ordinary Q em = — 1/3 quarks mix with the d! . 
Because the d' L has a different I^l from d^ Sl and b^, FCNC's appear at tree level 
in the left-handed sector. In particular, a Zbs coupling can be generated: 

£ Z FCNC=- 9 a U*S'fP L bZr + h.C. (3.1) 
/ COS fw 
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This coupling leads to a NP contribution to B s -B s mixing at tree level. 



The SM contribution to M[ 2 is 



Mtf M = ^ 2 M 2 w \ 2 t VBB B JlM B E(xt) 



(3.2) 



where X t = V t bV t * s , x t = ntf/M^, Vb — 0.551 is the QCD correction, and we take the 

1 /2 

hadronic parameter f Bs B B = 295 ± 36 MeV. The loop function E(x t ) is given by 



E(x t ) 



— 4x t + llx 2 — x\ 



+ 



3a;? In x t 



4(1 -xt) 2 2(1 -x t ) 3 

The mass difference AM S in the ZFCNC model is given by |28[ 

2 



(3.3) 



AM S = AM S 



SM 



i ■ I — ) -M^ 



u. 



sb 



X t 



Xt 



(3.4) 



where 



a = 4 



C(s t ) 



2^ 



^E(x t ) ' GpMyyElxt) ' 

The loop functions E(x t ) and C(xt) are given by p8[ 

-4xt + 1 lx 2 — xf 3xj In 



(3.5) 



C(x t ) = % 



4(1 -x t ) 2 2(1 -xt) 

4 — xt 3xt In xt 
4 [l-x t (l-xt) 2 



(3.6) 



The term in Eq. (|3.4| ) proportional to a is obtained from a diagram with both SM 
and NP Z vertices; that proportional to b corresponds to the diagram with two NP 
Z vertices. 

This coupling will also lead to contributions to T\ 2 due to the generation of the 
Z- mediated operator b — »• sf f (f = t, c). The amplitude is 



^ U sb s lfi (l - j 5 )bfY(I 3 - Q sin 2 9 W )(1 T 7s)/ ■ 



(3.7) 



There are diagrams with two NP vertices, and with one SM and one NP vertex. The 
dominant contribution is due to SM-NP interference, yielding 



r 



s,Z 
12 



G 2 F mlX c U s b r — 

— Vl - 4x c x 



7T 



1 



'1 - X r 



sm 2 9 w (l + 2x c ) } (0 LL ) 



+ 



sin 2 ^ {l + 2x c )(0 RR ) 



(3.8) 
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3.2 Z'-mediated FCNC's 



In the model with Z'-mediated FCNC's (Z'FCNC), the gauge group contains an 
additional U(l)', which leads to a Z' | 2"9|| . Within some string-construction or GUT 
models such as Eq, it is possible to have family non-universal Z' couplings. In the 
physical basis, FCNC's generally appear at tree level in both the left-handed (LH) 
and right-handed (RH) sectors. In particular, the interaction Lagrangian can contain 
Z'bs couplings: 



L z ' 

*-FCNC 



cos 9 



w 



S-fP L B%,b+(L<rHi) 



Z' + h.c. 



(3.9) 



which lead to a tree-level contribution to B s -B s mixing. (Note: the U(l)' coupling 
constant is g'. However, the above Lagrangian is written in terms of the SM coupling; 
the ratio of the two couplings is absorbed into the B s bS. This makes the comparison 
of NP and SM effects more straightforward.) 

The effective Hamiltonian for M[ 2 in the Z'FCNC model at the scale mj is 
calculated from the tree-level diagrams J3(| [31] : 

n cS = ^r^r V° /AS (B sb 

4 



4GV m| 
V2 m|, 



L " 2 LL (m b ) + rfl 2 \BlfO RR {m h ) + 2r ] 3 / 23 B^ b O LR (m b ) 



+ 3 (V 3/23 



(3.10) 



where rj = a s (M Z ')/a s {mb) and 

LL = 8^P L b8-fP L b. 
Orr = S llx P R b snfP R b 

LR = s^P L bs^P R b 
LR = sP L b sP R b . 



(3.11) 



In addition to the SM operator Oll, the operators Orr and Olr are generated at 
the scale Mz>- The operator Olr is generated through renormalization down to the 
scale m&. 



The matrix element for B s -B s mixing is given by 



M 



s,Z> 
12 



{B s \H c r\B s ) . 



The hadronic matrix elements are again calculated using the vacuum insertion ap 
proximation fl6|| . (B s \0 ll\B s ) is given in Eq. ( p.2|) , (Orr) 



(3.12) 
>n { 

(B s \Orr\B s ) 



(B s \O ll \B s ), and 
(Olr) 
(0 LR ) 



(B s \0 L r\B s ) 
(B s \Olr\B s ) 



1 
4 
1 

24 " 



1 

" 6 
1 

4 



m b + m t 



m B s f_ 



Tfifj + m s 



m Bs fl 



(3.13) 
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The Z' contribution to Y\ 2 through r pairs is obtained from the formulae given 



in Refs. ||, |33 



12 



AG 2 F ml ( m 2 7 \ 2 n — 



37T V m %' 

{B LL f + {B RR f + {B LR f + {B RL f][{\ - x T )(0 LL ) + (1 + 2x T )(0 RR ) 
+ 2 \B LL B RL + B RR B LR }\(1 - x t )(O lr ) + (1 + 2x T )(d LR ) 



+ q{b ll B lr + B RR B RL }x T (0 LL ) 
+ q{b ll B rr + B iK J B iiX }x r (0 LJi > 



(3.14) 



where B % i = B\ b B 3 TT = L,R), x T = rn^jm^. Comparing Eqs. ( |2.1| ) and ( |3.14j ), 
we can see that if B^ ~ (A c m|,/m|) 2 ~ 0(1), can be comparable to r^ SM - We 
will see that these Oil) couplings are still allowed by B(B S — > t + t~) and B(B — > 

X S T+ T-). 

Although the c-quark loop contribution with both couplings from NP cannot be 
large enough to compete with the SM contribution, there can be an interference term 
between the SM and the NP in this case. This term is given by 



12 



7r \ m z , 



AB^ b B L cc [{\ - x c )(0 LL ) + (1 + 2x c )(0 LR )} + QB^B r x c (Oll) + 12B R B R (0 LR ) 



+ AB R B L CC [{1 - x c )(0 LR ) + (1 + 2x c )(0 LR )} 



(3.15) 



Here again Eq. ( |3.15| ) can be comparable with r^ SM when B % sh B{ c ~ 0(1) 
L,R). 



4. Numerical Results 
4.1 ZFCNC 

In the numerical study of the ZFCNC model, we scan the allowed values of the 
couplings after imposing the constraints from |A S | = AM s /AMf M [Eq. ( p.4|) j and 
B(B — > X s fi + yr) 1 . The general expression for the branching ratio of B — > X s l + l~, 

1 There can be other constraints, such as those from the forward-backward asymmetry in B — > 
K* l + l~ |35|]. However, we do not include these in our analysis, since, as can be seen below in 
Fig. ||, the two constraints used are sufficient to rule out the ZFCNC model as an explanation of 
and a* v 
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where I = e, /i, r, is given by p4 

p(l-m s /m b ) 2 

B(B -> X s /+ /-) = B / [S 5Jlf (z) + B SM . VA {z) + By^)) cfe . (4.1) 



The expressions for £? , -Bsat, .B^m-vA; and .By a are given in Ref. p4| . We do not 
repeat them here. We note only that the integrand depends on the NP couplings 
R v ,Ra, R'v and i^. In the ZFCNC model, the NP couplings R' VA = 0, and the 
Rv,a are given by |36fl 



27T t/ s6 / 1 , „ . 2/1 \ „ 7T C/ Sl 



i?y = —-^ -- + 2sin^ w , R A = --^- (4.2) 
a K 6 y t * V 2 J aV tb V t * s 

We define g M as the sum of the 4- momenta of the /i + and \x~ in B — > X s ii + . 
The theoretical predictions for the branching ratio of this decay are not reliable 
over the whole q 2 region due to the presence of charm resonances at intermediate 
q 2 (7 GeV 2 < q 2 < 12 GeV 2 ). The predictions are relatively more robust for low 
q 2 (lGeV 2 < q 2 < 6 GeV 2 ) and high q 2 (14.4 GeV 2 < q 2 < m 2 ). However, the 
two regions have different sensitivities to the short- distance physics. The dominant 
contribution to B — > X s fi + fi~ in the low-g 2 region comes from a virtual photon, 
whereas the dominant contribution in the high-g 2 region comes from the Z and W. 
Since we are interested specifically in the Zbs coupling, we use the branching ratio 
in the high-g 2 region to constrain the ZFCNC parameter space. 




Figure 1: The 4>u ab -\U s b\ parameter space allowed by |A S | (magenta) and B(B — > 
X s [i + yr) (green). 

As described earlier, U s b denotes the Zbs coupling generated in the ZFCNC 
model [Eq. (|3.1| )1. The parameters of the model are therefore the magnitude and the 
phase of this coupling, \U s b\ and <fiu si = arg[/ s b. Fig. [I] shows the allowed region in 
(4>u sb , \U s b\) space due to the measurements of |A S | (magenta) and B[B — > X s /i + /i~) 
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(green). Using constraints only from | A s |, we see that the full range of <ftu sb is obtained 
only when \U s b\ < 0.001, while \U s b\ can be as large as 0.005 in the constrained regions 
of <pu Bb - However, if we include the constraints from B[B — > X s ji + the (4>u Bb , 
\U s b\) space is reduced considerably. From Fig. [I], one sees that B(B — > X s fi + 
limits the upper value of \U s b\ to be around 0.001. 

Using these allowed couplings, we can calculate the observables S^, afj and Ar s , 
and compare these values with the measurements. The la ranges of these observables 
are 



a sl 



[-0.91,-0.47] , 
[-17.7,-7.7] x 10~ 3 
±[0.084,0.208] ps" 1 . 



(4.3) 



In addition, we can examine correlations between the observables. 
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Figure 2: The a^-S^ correlation plot (left panel) and the AI^-S^ correlation plot (right 
panel). The horizontal lines in the right panel indicate the la experimental allowed range 
for Ar s . 

As a consequence of the strong constraint on the size and phase of U s b, we see 
in Fig. |2] that it is not possible to explain the present la ranges of and in the 
ZFCNC model. Still, this model can produce values for these observables which are 
closer to the experimental data. In particular, one can have ~ —3 x 10~ 4 , which 
is about 10 times larger in magnitude than the SM prediction of ~ 2 x 10~ 5 . And 

can be decreased to —0.08, as compared to the SM value of —0.04. 

4.2 Z'FCNC 

As discussed earlier, a large NP contribution to r^ NP , and hence to a s sl , may come 
from the NP b — > st + t~ and/or b — > sec operators. Since the two operators need 
not be related to one another, a large NP effect on b — > st + t~ but not on b — > sec, 
or vice versa, is conceivable. Thus, in analyzing the Z'FCNC model, we examine 
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whether it is possible to explain the present la measurement of a s sl if only one of the 
NP operators is present. 

4.2.1 NP b — y st + t~ operator only 

We first consider the case where the NP contribution to T S {2 P comes only from the 
b — > st + t~ operator. In our numerical study we scan the allowed values of the 
Z'FCNC couplings imposing the constraints from |A S |, B(B — > X s t + t~) < 5% 
and B(B S — > t + t~) < 5%. |A S | can be obtained using Eqs. ( [3.10| ) and (|3.12|) , 
and B(B — > X s t + t~) is given by Eq. ( [4.1| ) with I = r. The branching ratio of 
B s r+ T~ is |33| 



G%a 2 m 5 R fl t b , ,„ / 4m 2 

b(b s -+ t+t-) = 6 *;r Bs \v tb v;f > 



m. 



— J-{C\Q + Ra — R'a, 



m 



B, 



(4.4) 



In the Z'FCNC model, the NP couplings are given by 

47r ro|/m|, L L R 

K V T/ rr* B sbK B TT + B rr) , 

a v tb v ts 
47r m|/m|, L L R 

^ = - ! f Z ^5f b (^ T -^ T ). (4.5) 
a V tb V ts 

These couplings are constrained by the above observables. However, the key point 
is that, because the bounds on the NP b — > st + t~ transition are still weak, the 
constraints are not severe. 

Fig. [3] shows the a^-S^ and ATg-S^ correlation plots for the case where only 
Bg b and B^ R couplings are present. We make the following observations: 

• can take any value between —1 and 1. 

• The Z'FCNC model gives values for a s sl within its la range. 

• One can simultaneously reproduce the present measurements of all three ob- 
servables, S^, a s sl and Ar s , within their la ranges. 

Obviously, one obtains the same conclusions if a nonzero Bf b is also present. 
4.2.2 NP b —> sec operator only 

We now consider the case where the NP contribution to T*2 NP comes only from the 
b — > sec operator. In our numerical study we scan the allowed values of the Z'FCNC 
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Figure 3: The a^-S^^ correlation plot (left panel) and the Ar^-S^ correlation plot (right 
panel), for the case where only Bri and B T r couplings are present. The horizontal and 
vertical lines indicate the la experimental allowed ranges for the observables. 



couplings imposing the constraints from |A S | and the indirect CP asymmetry in 
Bd — > J/tpKs- This latter constraint arises as follows. In the presence of NP in the 
decay b — > sec, the effective measured sin 2(3 in Bd — > J / ipKs is given by [ p3[ 



sm2^ meas = sin2/3 + 2|r|cos2/3sin^cos(5 a . (4.6) 



Here, the NP parameters are defined similarly to those in Eq. ( |2.10| ): |r| is the ratio 



of magnitudes of the NP and SM b — > sec amplitudes, ip is the NP weak phase, and 
5 a is the SM strong phase. The true value of sin 2(3 is taken from the fit to the sides 
of the unitarity triangle: sin 2(3 = 0.731 ±0.038, while the experimental measurement 



gives sm2(3 meas = 0.668 ± 0.028 [p7|| . Using these values, |r| is estimated to be [p3 



\r\ = (4.6 ±3.5)%. (4.7) 

Thus, at la, \r\ < 8.1% is permitted. This upper bound on \r\, along with the 
measurement of |A S |, constrain the NP couplings. 

Fig. |] shows the a s sl -S^<p and ATg-S^ correlation plots for the case where only 
B^ b and B^ R couplings are present. The results are the same as for the case where 
a NP b st + t~ operator was added: 

• S 1 ^ can take any value between —1 and 1. 

• The Z'FCNC model gives values for a s sl within its la range. 

• One can simultaneously reproduce the present measurements of all three ob- 
servables, S^, a s sl and Ar s , within their la ranges. 

We therefore see that the Z'FCNC model can reproduce the present measure- 
ments of 5^,0, Ar s and within their la ranges if the NP contribution comes only 
from either the b — > st + t~ or b — > sec operators. 
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Figure 4: The a^-S^^ correlation plot (left panel) and the Ar^-S^ correlation plot (right 
panel), for the case where only Bri and B C c couplings are present. The horizontal and 
vertical lines indicate the la experimental allowed ranges for the observables. 



5. Conclusions 

Recently, the D0 Collaboration measured the CP-violating like-sign dimuon charge 
asymmetry in neutral B decays and found a 3.2<r difference from the prediction of the 
standard model (SM). Combining the D0 result with other measurements, we obtain 
the semileptonic charge asymmetry for B s mesons, (a^ave = —(12.7 ± 5.0) x 10" 3 , 
which is still ~ 2.5a away from the SM prediction. 

A non-SM can be explained theoretically by a new-physics (NP) contribution 
to B s -B s mixing, specifically to M[ 2 and r* 2 , the dispersive and absorptive parts of 
the mixing. It is usually assumed that rf 2 is dominated by the SM tree-level b — > sec 
coupling, and that it is not affected by the NP. Instead, NP is considered only in 
Mf 2 . However, in this case, the measured value of (a^ave within its la range cannot 
be explained. Thus, NP in is necessary. 

There are some NP models (leptoquarks, R-parity-violating SUSY) whose con- 
tributions to rf 2 can compete with the SM. In this paper, we examine models with 
flavor-changing neutral gauge bosons {Z or Z') to see if they can explain (a^ave- 
The models must also reproduce the measured values of the indirect CP asymmetry 

in B s — > J/ip(j) and AT S . 

In our models we assume that the main contribution to T* 2 comes from diagrams 
involving r and/or c loops. The point is that the current constraints on NP b — > sff 
(f = t, c) transitions are quite weak. We find that the model with flavor-changing 
Z couplings (ZFCNC) cannot reproduce the present measured values of and 
within their la ranges. Still, the ZFCNC model can lead to modified values for 
and which are closer to the experimental data than the SM predictions. 

On the other hand, in the model with flavor- changing Z' couplings (Z'FCNC), 
even after imposing the current experimental constraints, the present measurements 
of S^, a s sl and Ar s , can all be simultaneously reproduced within their la ranges. 
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Indeed, the full range of is allowed. We therefore see that the Z'FCNC model is 
another type of NP which can explain the measured value of the D0 dimuon charge 
asymmetry in the B system. 
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